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Swirling-Sweepers: 
Constant- Volume Modeling 


Alexis Angelidis* Marie-Paule Canit 
University of Otago Laboratoire GRAVIR* 
Abstract 


Swirling-sweepers is anew method for modeling shapes 
while preserving volume. The artist describes a deforma- 
tion by dragging a point along a path. The method is in- 
dependent of the geometric representation of the shape. 
It preserves volume and avoids self-intersections, both lo- 
cal and global. It is capable of unlimited stretching and the 
deformation can be constrained to affect only a part of the 
model. 


We argue that all of these properties are necessary for 
interactive modeling if the user is to have the impression 
that he or she 1s shaping a real material. Our method is 
the first to implement all five. 


1. Introduction 


In a virtual modeling context, there is no material: no 
wax, Clay, wood or marble. A challenge for computer 
graphics is to provide a virtual tool that convinces the 
artist that there is material. To perfect this illusion, the 
shape must behave in accordance with a suitable mod- 
eling metaphor. 


Volume is one of the most important factors influenc- 
ing the manner in which an artist models with real ma- 
terials. A virtual tool preserving volume is needed to 
help the artist believe he is interacting with material. 
Also, modeling by preserving the available amount of 
material will produce a shape with style, that other vir- 
tual modeling methods can only achieve with more ef- 
fort. 
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1.1. Previous volume-control models 


Volume preservation has been recognized for a long 
time in animation, as a desirable property for the ani- 
mation of believable animal and human characters [13]. 
[9] use constrained optimization methods for objects 
discretized into lattices. [5] use controllers for main- 
taining the implicit surface that coats a set of particles 
to a constant volume during deformation. [7] achieve 
incompressibility in water simulation by maintaining 
a divergence free velocity field, thanks to the Poisson 
equation. [12] rely on finite volume methods to sim- 
ulate quasi-incompressible materials such as muscular 
tissue. 





Volume preservation has also been considered as a very 
useful constraint for the intuitive modeling of shapes. 
[10] propose an optimization method to adjust the con- 
trol points of the popular free form deformations (FFD) 
[11], but it works only for tensor-solids. [8] also adjust 
FFD control points, but their method does not allow 
local editing. [3] propose a volume preserving space de- 
formation based on a model called DOGME. The defor- 
mation does not have a local support, and requires the 
computation of the shape’s volume. [4] preserve only a 
volume between the surface and a base surface. [6] in- 
troduce mass-preserving local and global deformations 
for shapes represented by a mass-density field sampled 
in a grid. 





The limitation of existing methods is either that they 
only apply to a specific type of geometric representa- 
tion, or they only apply to shapes whose volume can 
be computed. 


1.2. Overview 


This paper presents swirling-sweepers, a new method 
dedicated to modeling shapes while preserving the 
shape’s volume. Our technique belongs to space de- 
formations, and is therefore applicable to a wide range 


of geometric representations, including all of the pop- 
ular parametric surfaces. 


It is the first method that preserves volume, has a lo- 
cal support, prevents local and global self-intersection 
of the surface and does not require any volume com- 
putation. Most importantly, using the method is sim- 
ple: the artist only has to provide the trajectory of a 
point, for instance with a mouse. 


In Section 2 we summarize the principle of the space 
deformations called sweepers. Then we present in Sec- 
tion 3 our new method for modeling by constant vol- 
ume deformation. 


2. Principle of Sweepers 


We briefly review the elements required for understand- 
ing the space deformations called sweepers [2]. 


Space deformation provides a formalism to specify any 
modeling operation by successively deforming the space 
in which an initial shape, S(to), is embedded. A de- 
formed shape is given by the modeling equation’: 
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where ftot : R? — R° are n space deformations, de- 
forming a point, p, of shape S(t;) into a point of shape 
S(ti+1). 


Informally, a sweeper is a geometric tool together with 
a motion path. This tool defines an influence function. 
The basic idea is that the tool is placed somewhere in 
the region of a shape to be deformed and moved along 
the path. The motion drags a part of space defined by 
the influence function, in a manner that prevents the 
shape from self-intersecting. 


More formally, the deformation is defined by a scalar 
function, ¢; : RÌ + [0,1], that varies over time, t. This 
field is defined by composing the distance to the tool, 
d+, with an influence function, u 


De = ode (2) 


Any smooth decreasing function of finite support can 
be used for u. We use a C° continuous piecewise poly- 
nomial, in which À defines the radius of the influence 
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A deformation is defined by transforming the tool’s po- 
sition, size and orientation, given by the matrix M, 
into the next configuration, given by the matrix M;,,,. 
Let us denote M; = M;,,, M; the transformation ma- 
trix from the previous to the new configuration. A naive 
deformation of a point p with a single tool would be 


be (p) = (de, (p) © M;) p (4) 


where the matrix operator © is defined asa © M = 
exp(alog M). In Section 3.3, we give the closed-form 
for computing exp and log. We refer the reader to [1] 
for a more detailed overview. Loosely speaking, the op- 
eration © is the equivalent of multiplying a matrix by a 
scalar. It raises a transformation to a non-rational ex- 
ponent. However, Equation 4 does not prevent the sur- 
face from self-intersecting. By decomposing the defor- 
mation into s sub-functions [2], self-intersections are 
avoided 


s—1 
Tiaa un Voa (p) (5) 
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For a tool in a bounding volume Vp, the number of 
steps is 


s = max(1, | — min(ôux/0d)max og M;p||]) (6) 
POVb 


3. Constant Volume deformation 


We introduce here swirling-sweepers as a method 
for deforming shapes while preserving their vol- 
ume. Swirling-sweepers are a particular case of sweep- 
ers that use only point tools. 


3.1. A basic deformation 


We define a particular case of sweeper, a swirl, by us- 
ing a point tool, c, together with a rotation of angle 0 
around an axis U (see Figure 1). A scalar function, ¢, 
and a deformation are defined as before (see Equations 
2 and 4). Informally, a swirl twists space locally around 
axis v without compression or dilation. We prove in Ap- 
pendix A that a swirl preserves volume. 





3.2. Combining for complexity 


Many deformations of the above kind can be naively 
combined to create a more complex deformation 


f(p) = [Bow © 1) p (7) 


i=0 


where @ is a commutative addition of transformations 
defined as M @ N = exp(log M + log N) [1]. We pro- 
vide a convenient way for the artist to input n rota- 
tions, by specification of a single translation t. Let us 
consider n points, c;, on the circle of center h, and ra- 
dius r lying in a plane perpendicular to t. To these 
points correspond n consistently-oriented unit tangent 
vectors v; (see Figure 2). Each pair, (ci, 0i), together 
with an angle, 6;, define a rotation. Along with radii of 
influence A; = 2r, we can define n swirls. The radius of 
the circle, r, is left to the user to choose. The follow- 
ing value for 6; will transform h exactly into h +t (see 
Appendix B). 








jt |] 
Ô: = 8 
i nr (8) 
With this information, the deformation of Equation 7 
is now a tool capable of transforming a point into a de- 
sired target. We show in Figure 2 the effect of the tool 
for different values of n; in practice, we use 8 swirls. 





Preserving coherency and volume If the magnitude of 
the input vector tis too large, the deformation of Equa- 
tion 7 will produce a self-intersecting surface, and will 
not preserve volume. The reason for self-intersection is 
explained with details in [2]. The volume is not pre- 
served because the blending operator, @, blends the 
transformation matrices, and not the deformations. To 
correct this, it is necessary to subdivide € into smaller 
vectors. ‘The number of steps must be proportional to 
the speed and inversely proportional to the size of the 
tool. We use 








s = max(1, [4|lé ||/r]) (9) 


As the circle sweeps space, it defines a cylinder. Thus 
the swirling-sweeper is made of n - s basic deforma- 
tions. Figure 3 illustrates this decomposition applied 
to a shape. 


3.3. Swirling-sweepers algorithm 





We summarize here the swirling-sweepers algorithm: 
Input point, h, translation, t, and radius, r 
Compute the number of required steps, s 
Compute the angle of each step, 0; = 2iléll 

for each step 7 from 0 to s — 1 do 
for each point p in the tool’s bounding box do 

M=0 

for each swirl à from 0 to n — 1 do 
M+= par (|p — cij||) log M, 

end for 

p = (exp M)p 


end for 
end for 


h 


The point c;; denotes the center of the in swirl of the 


jth ring of swirls. For efficiency, a table of the basic- 
swirl centers, c;;, and a table of the rotation matri- 
ces, log M; ;, are precomputed. We have a closed-form 
for the logarithm of the involved matrix, given in equa- 
tions 10 and 11, saving an otherwise expensive numer- 
ical approximation. 


no == 0, 
m= Ci j X n (10) 
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Note that for the sake of efficiency, we handle these ma- 
trices as mere pairs of vectors, (fi, m). Once M is com- 
puted, we use a closed-form for computing exp M. Since 
the matrix M is a weighted sum of matrices log M; ;, 
the matrix M is of the form of Equation 11, and can 
be represented with a pair, (ñm, ™m). If ñm = 0, then 
exp M is a translation of vector mpm. Else, if the dot 
product Mu - Nu = 0, then exp M is a rotation of cen- 
ter €, angle 0 axis v, as given by Equation 12. 
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e = 
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0 —=[ñwl 
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Finally, in the remaining cases, we denote l = [nl 


and we use Equation 13. See Appendix C for efficiency. 
expM=I1+M+>9M?+58™M? (13) 


Symmetrical objects can be easily modeled by intro- 
ducing a plane of symmetry about which the tool is re- 
flected (see Figure 5). 


4. Results 


We have implemented swirling-sweepers in C++ using 
OpenGL®, on a Pentium® 2400Mhz with 1GB of RAM. 
This implementation works in real-time. The compu- 
tational time is a function of the magnitude of the in- 
put vector, because this determines the number of sub- 
steps. Small vectors will produce extremely fast defor- 
mations. In order to preserve the sampling of the de- 











formed surface, we use the mesh update algorithm pro- 
posed in [2], adapted for sweeping space deformations. 
Results are shown in Figure 5. 


Limitations In our implementation, the tool must be 
of significant size compared to the density of the mesh. 
In Figure 5, we compare the shapes’ volume with unit 
spheres on the right. The shapes volumes are respec- 
tively 101.422%, 99.993%, 101.158% and 103.633% of 
the initial sphere. This error is the result of accumu- 
lating smaller errors from each deformation. For in- 
stance 80 swirling-sweepers have been used to model 
the alien. The small errors are due to the finite num- 
ber of steps, and to our choice of shape representation. 


5. Conclusions and future work 





We have presented swirling-sweepers, a new volume- 
preserving space deformation that uses the sweepers 
formulation: a combination of matrices raised to pow- 
ers of scalar functions. Combined with the original 
sweepers, the volume of a shape can be increased, pre- 
served or decreased. We believe there are many more 
useful sweeper operations yet to be discovered, for in- 
stance more complex volume-preserving tools, topol- 
ogy changing tools, surface area preserving tools, or 
surface smoothing tools. We also believe our technique 
is adaptable to volume preserving animation. 
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A. Constant volume basic swirl 


Let us prove that a basic swirl preserves the volume. 
We will show that the determinant of the Jacobian, 
det J, is equal to 1 everywhere. A rotation of angle 20 
around n and its powers in @ can be modeled with a 
quaternion: 


q? = (cos(#8), sin(p0)ñ) (14) 
(x, y,z)* is 
F(p) = apg? (15) 


We can assume without loss of generality that the ro- 
tation is centered at the origin. To express the Jaco- 
bian, we need the three partial derivatives of f. Let us 
denote ¥, ÿ and Z the vectors (1,0,0)7, (0,1,0)' and 


Hence the swirl deformation of a point p = 





(0,0,1)!'. The first partial derivative of f is 
df(p) _ 8q? bq? 
ag? 09 + Gf 1 
eae 0 a (16) 
The reader can verify the quaternion equality: 
ôq? ob 5 
eE — (sin(0), — cos(ô)n) (17) 


Using 17, the leftmost term of the right side of equa- 
tion 16 is a quaternion: 


Be ba" ge 
= + sin(b0), cos(b@)n) * p * q? 
— 922 (— cos(p0)n - p, — sin(0)p + cos(p0)n x p) * q? 
= 932 (— ñ -p,cos(290)(n x p) — sin (290) (ñ x p) x ñ) 


The reader can verify similarly that the rightmost term 
of equation 16 is also a quaternion: 


5a” pq? = OSE (R - p, cos(2¢0) (7 x p) — sin(260) (7 x p) x 7) 


The partial derivative of f in x is a vector: 
Sf) — 9952 (cos(20) (7 x p) — sin(2G0)(# x p) x R) + g?Xq? 


Let us introduce g, and ü for the sake of simplicity: 


qz = q? xq? 
u = 20(cos(290)(n x p) — sin(2¢0)(7% x p) x 7) 


The partial derivative in x shortens to: 


df(p) _ ô$ 

———© = U +, 18 
The two other partial derivatives of f are obtained by 
substituting x for y or z. Since a rotation preserves 


lengths and angles, we can write: 


> 


Ux = Ty X q: 


Let us develop the determinant of the Jacobian: 


= “OF COR NOR 


= 76 ô ô 
= 1+. (gr + Lay + a) 

We can assume without loss of generality that n = X. 

This provides expressions for the rotated canonic set: 


> 


qz = A 
Ty cos(290)ÿ + sin(2¢0)Z 
TA cos(290)z — sin(2¢0)¥ 


This assumption also provides a simple expression for 
the double cross product: 


(ix p) x n = yy + zZ 


We will now use the fact that the tool is spherical. We 
model the field function @ as a function of the distance 
to the origin, d(p). The field can be partially derived: 


6O(d(p)) _ db dd(p) ô x 





ôx dd dx  ddd(p) 
With this, the determinant of the Jacobian becomes: 
det J 
= Ata 


(cos()X x p — sin()(y¥ + 22)): 
(aX + (y cos() — zsin())¥ + (ysin() + z cos())Z) 
1 


Thus det J is equal to 1 everywhere. Therefore the de- 
formation stretches space with no expansion nor com- 
pression. 





B. Swirl angle 


The image of a point p in the center of a circle of swirls 
is given by Equation 7. Since the point is at the center, 
one can substitute ġ; for 1/2 

f(p) = Pio (3 © Mi)p (19) 
The speed of this deformation at p is given by the log- 
arithm 


B= > à log(Mi)p (20) 


Since M; is a rotation matrix, this simplifies (see Equa- 
tion 11). 





C25 =e) (21) 
By taking the norm: 
lel = 3 Dio ll — cill (22) 
Since the centers are equidistant to p 
ol = Sr (23) 





Therefore the angle is 
g = 2151 


nr 


(24) 





C. Exponential 


Applying the exponential of the matrix to a point does 
not require to compute the exponential of the matrix 
explicitely. Let us define the matrix M with a pair of 
vectors, (n, m). 





exp(M)p = p+(m+nx p)b+ ("$> — pla 
+ A(T x pja + (mx Mm)(1 — b))z 
where 1 = ||7|| (25) 
a = 1-—cos(l) 
—  sin(l) 
pa n 
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Figure 1. The effect on a sphere of a swirl cen- 
tered at c, with a rotation angle 0 around v. The 
two shapes have the same volume. 





Figure 2. By arranging n basic swirls in a circle, 
a more complex deformation is achieved. In the 
rightmost image: with 8 swirls, there are no visi- 
ble artifacts due to the discrete number of swirls. 
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Figure 3. A volume preserving deformation is ob- 
tained by decomposing a translation into circles 
of swirls. 3 steps have been used for this illus- 
tration. As the artist pulls the surface, the shape 
gets thinner. The selected point’s transforma- 
tion is precisely controlled. 





Figure 4. When pushed or pulled, a sphere will 
inflate or deflate elsewhere. 











Figure 5. Examples of models “sculpted” with 
swirling-sweepers. The mouse, the goblin, the 
alien and the tree have respectively 27607, 
25509, 40495 and 38420 vertices. These objects 
were modeled in less than 30 min by one of the 
authors. Eyeballs have been added. 


